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Ahtnct 

A  simple  epidemic  process  in  which  the  number  of  individuels  who  can 
become  infected  at  any  point  in  time  is  itself  a  random  variable  is  described.  The 
discrete  asymptotic  behaviour  of  such  a  process  is  discussed.  In  particular,  the 
associated  marfinal  distribution  of  the  limiting  process  is  considered. 

SIMPLE  BATCH  tnonc  ntOCESS:  UMmNO  DISTIUBLmONS;  MARGINAL  DISTRIBU¬ 
TION 


1.  Introduction 

A  simple  epidemic  process  in  a  closed  population  consists  of  susceptible 
individuals  who  may  become  infective  individuals.  Suppose  initially  there  are  N 
susceptibles  and  a  infectives.  Recently,  Biltard,  Lacayo  and  Langbcrg  (1979) 
showed  how  this  process  may  be  described  in  terms  of  N  independent 
exponential  random  variables  representing  the  times  between  any  two  successive 
occurrences  of  a  new  infective. 

In  this  work,  the  following  extension  of  the  simple  epidemic  process  is 
considered.  When  occurrences  of  infectives  in  the  simple  epidemic  model  are 
viewed  as  realizations  of  a  jump  process,  each  jump  is  of  fixed  size  one. 
However,  more  generally,  the  size  of  the  jumps  may  be  a  random  variable.  For 
descriptive  reasons,  we  call  this  extension  a  simple  batch  epidemic  process. 

Our  purpose  is  to  investigate  the  asymptotic  behaviour  of  a  sequence  of  simple 
batch  epidemic  processes  when  the  population  size  N  goes  to  infinity.  Before 
obtaining  the  main  results,  we  elaborate  and  introduce  some  notation. 
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2.  A  simple  batch  epidemic  model 

Let  the  size  of  the  jumps  be  represented  by  the  sequence  of  random  variables 
Zh  i  =  1,2,  •••,  with  Z,  taking  values  {1,2,  •••}.  That  is,  Z  is  the  number  of 
individuals  becoming  infected  at  the  it h  occurrence  of  infection.  Then,  the 
number  of  jumps  that  has  occurred  until  at  least  k  of  the  susceptible  individuals 
are  infected,  is  R(k)  given  by 

(1)  J?(k)“  min  jrj 2  Z,  S  fcj,  fc  =  1, 

Clearly,  1  £  R(k)£  k.  If  D,  is  the  total  number  of  new  individuals  who  have 
become  infected  just  before  the  ith  jump,  it  follows  that 


(2) 


0.  i  =  1, 

Sz;  i  =2,*--,R(N), 

i- 1 


where  initially  there  are  a  infectives.  Thus  (2)  implies  that  after  the  ith  jump 
{OS  i  £  R(N)~  1]  there  are  (D,„  +  a)  infectives,  and  after  the  f?(lM)th  jump, 
there  are  (N  +  a)  infectives.  For  convenience,  we  take  a  =  1. 

We  recall  that  in  the  simple  epidemic  process  (where  Z,**  1  always),  the 
interinfection  time  T,  between  the  (i  -  I)fh  and  ith  occurrence  of  a  new  infective 
is  exponentially  distributed  with  parameter  proportional  to  i  (N  +  1  -  i).  How¬ 
ever,  in  the  present  case,  since  the  jump  sizes  are  themselves  random  variables, 
the  conditional  mean  (interinfection)  time  is  also  a  random  variable.  Therefore, 
let  us  define  the  realizations  of  the  random  variable  p(N,  D(),  <■*!,••  • ,R(N ), 
as  the  mean  time  between  the  occurrence  of  the  (i  -  l)th  and  ith  jump,  that  is, 
E(Tt). 

We  may  now  define  the  interinfection  time  between  the  (i  -  l)th  and  ith  jump 
as  T,  » (i(N,D,)Uh  i «  1,  •  •  •,  R(N),  where  the  U„  i « 1,  •  •  •,  R(N),  are  inde¬ 
pendent  exponentially  distributed  random  variables  with  mean  one,  and  where 
the  If,  are  independent  of  the  Z*  /  «  1,2,  ■  •  It  follows  immediately  that  the 
total  time  until  there  are  at  least  k  new  infectives,  &,  may  be  written  as 

(3)  &  -  S*  m  W  D.)t/„  k  «  1,  •  •,  R(N). 

1*1 


We  note  that  the  particular  case  S«,w, ■  SN  is  simply  the  duration  time  of  the 
simple  batch  epidemic  process. 

Finally,  we  denote  the  total  number  of  new  infectives  present  at  time  t  by 
XN{t),  r  >0,  with  realizations  0,  •  •  •,  N.  Then,  for  every  set  of  positive  real 
numbers  th  /  ■!,•••,)»  and  a  corresponding  set  of  positive  integers  k„  j  = 
!,*•*,  m,  we  have  the  relationship 
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(4)  P<x*(4)  e**/  *  l,  ■••,«}.=  J»{&,  s  u  - 1.  *. «}. 

Hence,  it  follows  that  the  simple  batch  epidemic  process  is  uniquely  determined* 
by  the  p(N,  D,).  i  =  1,  •  •  *,  R(N),  and  the  Z*  i  -  1,2, •  • 

In  Section  3,  we  present  a  condition  which  is  sufficient  to  ensure  the 
convergence  in  law  of  a  sequence  of  such  processes  to  a  discrete  stochastic 
process.  In  Sections  4  and  S,  some  specific  marginal  distributions  of  X(t)  are 
investigated. 

3.  Convergence  in  law  of  a  sequence  of  simple  batch  epidemics 

Let  X(t),  l  > 0,  be  a  stochastic  counting  process.  Let  there  be  associated  with 
X(t)  the  random  variables  /t(Dt)  where 

0,  *  =  1. 

(5)  D,= 

«  =  2,3,-  •, 

with  the  Zf  being  the  jump  random  variables  defined  in  the  previous  section. 

Let  X„(t),  I  >0,  N  =  1,2,  •  ■  be  a  sequence  of  simple  batch  epidemics  as 
defined  in  the  previous  section.  On  the  array  of  functions  n(N,D,), 
i  *  1,  •  •  -,  R(N),  N  ■  1,2,  •  •  -,  we  impose  the  following  condition.  For  each 
value  of  i. 


(6)  iimM(N,0)  =  M(A)a.s. 


where  ft  (Dt),  i  -  1, 2,  •  •  -,  is  a  finite  positive  function.  One  such  family  of  double 
arrays  of  mean  interinfection  times  satisfying  (6)  is 

(7)  /t(N,  Du)  =  cNm(N  -  D,y(D,  +  1)"*,  i  =  1,---,R(N),N  =  1,2,—, 

where  a,  0,  and  c  are  positive  quantities.  This  family  may  be  viewed  as  a 
generalisation  of  that  suggested  by  Severo  (1969).  Another  particular  case  is  the 
model  of  McNeil  (1972). 

Let  the  waiting  times  between  the  (i  -  l)th  and  ith  jump  be  T\"  ft(D,)Uh 
i  =  1,2,-  •*,  where  as  before  the  U,  are  independent  exponentially  distributed 
random  variables  with  mean  one,  and  are  independent  of  the  Z„  j  =  1,2,---. 
Finally,  let  us  define 

(8)  M(Ot)Uh  *  =  1,2,  -, 

1-1 

where  R  (*)  is  as  defined  in  (1).  Then,  for  any  arbitrary  but  finite  set  of  positive 
real  numbers  t„  j  =  1,  •  •  -,  m,  and  positive  integers  fc„  /'  =  1,  •  •  -,  m, 
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(9)  P{X(t,)2  k,  j  -  1,  •  •  m)  -  P[Si,  S  W  -  1,  •  •  •. «}. 

The  following  theorem  may  now  be  proved. 

Theorem.  Let  X*(f),  t>  0,  N  »  1,2,  •  •  be  a  sequence  of  simple  batch 
epidemics  as  described  by  the  quantities  m  W  D, ),  i  •  1,  •  •  %  R(N),  and  Z* 
i  1,2, When  the  condition  (6)  holds,  the  joint  and  the  marginal  state 
probabilities  of  the  XM(t)  process  converge  to  the  corresponding  state  prob¬ 
abilities  of  the  counting  process  X(t),  t>  0. 

Am/.  We  prove  the  result  for  the  joint  state  probabilities.  The  result  for  the 
marginal  probabilities  is  proved  analogously.  It  is  sufficient  to  show  that  for 
every  arbitrary  but  finite  set  of  positive  real  numbers  /„/'=*  1,  •  •  -,m,  and  positive 
integers  k„  >  -  1,«  •  -,m, 

(10)  Km  P{XN(0  St  kf,i  "  1,  •••,»»}»  P(S;  S  4,/  -  1,  •  *  m>. 

To  prove  (10),  we  recall  that  R(k,)  is  independent  of  N  and  is  bounded  above  by 
kh  The  condition  (6)  implies  that 


Thus, 

S*  -  2  M (N,  D,)17,  2  m(A)1/,  -  Si, 

l-l  «»l 

tor  all  /  «  1,  •  •  •,  m.  Therefore, 

(11)  KmF(&,a4,/*  1,  •••,»»}  «P{S;,  a  4,/-  1,  •••,!»}. 

Hence,  combining  (4)  and  (11),  the  result  (10)  follows. 

We  note  that  the  exponentiality  assumption  of  the  interinfection  times  as  well 
as  the  independence  between  the  U,  and  the  Z%,  i  - 1,2,  •  *  *,  were  not  necessary 
to  prove  the  theorem.  Thus,  this  theorem  may  be  applied  to  more  general 
situations. 

4.  Distribution  of  the  limit  tag  epidemic  process 

In  this  section,  a  study  of  the  marginal  distribution  of  X(t)  is  made.  In 
epidemic  theory,  this  corresponds  to  the  distribution  of  the  number  of  new 
infectives  for  a  large  population.  It  follows  that 
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P{X(t)a  k}  -  p{"£  md,)U,  s  r) 

<12)  -  ij  pjs  m(D,)W  S  l|R(fc)-  v}p{R(k)-  v}. 

Cleariy,  we  may  write 

P{R(k)  -  *}  -  p(£  Z  <  k,  2i  £  k) 
and  if  the  Z  are  independent  random  variables 

(13)  P{R(k)-*}“  £  p(§Z“«)p(Z.iSk-n.). 

Therefore,  when  the  2  are  identically  distributed,  this  probability  (13)  can  be 
easily  evaluated. 

The  conditional  probability  term  of  (12)  may  be  written  as 
p{|;M(D,)f/(S<|R(k)-vJ 

(14)  ,  .  , 

-  £  m (A)«  S  l| R(k) -  v, Z  »  rj  P(Z  -  x), 

where  x  «  (z,t  •  •  •,  x.)  is  a  particular  realization  of  a  set  of  jump  statistics 
Z  »  (Z,,  •  •  •,  Z.)  and  where 

"*  1**4 

If  the  Z  are  independent  identically  distributed  random  variables, 

(15)  P(Z«x)-]Qp(Z-x,). 

To  evaluate  the  remaining  term  in  (14),  it  is  noted  that  for  given  Jt(k )  and  Z,  the 
actual  jump  sizes  at  each  stage  ate  known.  Hence,  the  distribution  of  p  (D,)U,  is 
known  for  each!  "  1,  •  •*,»>.  If  the  parameters  p( A),  i  ■  l,*",*  are  distinct,  an 
application  of  Billard,  Lacayo  and  Langberg  (1979),  Theorem  1  yields  the 
desired  solution;  whilst  if  the  parameters  p(A),  i m  1,  •  •  v,  are  not  distinct  the 
solution  follows  from  an  application  of  Billard,  Lacayo  and  Langberg  (1990), 
Theorem  2.  Thus,  we  obtain  the  distribution  function  from  (12)  and  hence  the 
probability  distribution  of  the  number  of  infectives. 

Although  we  have  restricted  our  attention  here  to  the  marginal  distribution  of 
X(t),  we  note  that  an  explicit  form  for  the  finite-dimensional  joint  distribution  of 
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the  process  {X(4),/  “  1.  •  •  -,m}  can  be  easily  studied  by  generalizing  (12),  (13) 
and  (14). 


5.  Examples 

There  are  some  particular  cases  worthy  of  special  attention.  Consider  first  the 
distribution  of  R(k).  Suppose  that  the  random  variables  (Z  - 1)  are  indepen¬ 
dent  and  follow  a  Poisson  distribution  with  parameter  A.  That  is, 

P{Z  =  z  + 1}  =  e-kk'/z\, z  - 0, 1, •  • 

Then,  from  (13)  it  can  be  shown  that 

P{fl(fc)-i'}«  2  2  e~k',km*’~¥*\v  -  iy"“*’*,/[r!  (m  —  + 1)!  J 

■t-a-1 

(16)  =  q(v),  say. 

Another  example  of  jump  processes  is  when  the  Z,  are  independent  random 
variables  taking  values  1  or  2  with  probability  p  or  q,  respectively.  This 
corresponds  to  the  case  in  which  individuals  can  become  infected  in  batches  of 
one  or  two  only.  Then, 

P{R(k)~  v)  -  p{§  Z,  -  k  - 1. Z.  =  1  or  2}  +  p{£  Zi  -  k  -  2, Z.  =  2} 

(17)  m(k-v)pti ,qk 'r(2v~k9- 

Let  us  now  consider  special  forms  for  the  quantity  One  particular  case 

is  when 

m(D,)=  Jimp  (N,D>)-$mNm(N -Dty  (ml)\ 
where  a,  fi  >0  and  where  m  *  E(Zi)>0.  That  is, 

(18)  M(D,)-(mi)*. 

The  equation  (12)  gives 

P{X(t)S  k} -  2  p{2  (miy'U,  5  t| R(k)  =  v}p{P(fc) -  v). 


and  for  0-1, 


(19) 


P(X(f)a  k)  *  i;  pf2  r'W  S  mil  •  P{*(k)-  v}, 

ir*l  ilmt  * 

=  Z(l-0'Wk)»»( 


by  an  application  of  BiUard,  Lacayo  and  Langberg  (1979),  Theorem  1. 
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In  some  populations,  if  N  is  sufficiently  large  relative  to  the  number  of 
infectives,  it  is  reasonable  to  expect  that  the  infection  rate  is  a  constant 
independent  of  the  number  of  infectives.  That  is, 

(20)  M(A)=limM(Ar,D,)  =  c. 

Then  it  follows  that 

pfct  l/iSf|R(k)=  v}p{R(k)-  v} 

w-l  l  i-J  i 

(21)  -  i  P(W.  S  c'r)P{R(*)  = 

a«l 

where  W.  has  a  gamma  distribution  with  parameters  1  and  v. 

Finally,  the  quantity  (k(D,)  may  still  be  a  random  variable  but  also  indepen¬ 
dent  of  the  number  of  infectives.  That  is, 

(22)  M(A)-D. 

Then, 

P{X(t)S  k}  -  2  P(Y  -  DW.  3  f)P{R(k)*=  v) 

•■1 

where  W.  follows  a  gamma  distribution  with  parameters  1  and  v.  For  a  given 
distribution  D,  this  probability  can  be  evaluated.  For  example,  suppose  D 
follows  a  Pareto  distribution 

/(d)«o0*d— *,  d>0,  0,0  >0. 

Then,  routine  but  lengthy  derivation  gives 

(23)  PCKsO-nvHHMW-Wl^  +  M/W, 

where 

r(*f€)mfc  e~muA~'du 

is  the  tabulated  incomplete  gamma  function. 

We  note  that  we  may  write 

where  N(t)  is  the  number  of  jumps  that  have  occurred  up  to  time  t.  If  p(D,)  «  c, 
N(l)  is  a  Poisson  process  with  rate  c~*  and  thus,  if  the  Z,  are  independent  and 
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identically  distributed,  X(t)  is  a  compound  Poisson  process.  II  the  ft(D,)  are 
independent  and  identically  distributed,  N(t)  is  a  renewal  process  and  thus,  if 
the  pain  (ZkM(A)]  are  independent  and  identically  distributed,  X(t)  is  a 
renewal  reward  process  (Ross  (1970),  p.  51),  that  is  a  special  case  of  a  cumulative 
process  (Cox  (1962),  p.  91). 


The  authors  wish  to  thank  the  referee  for  his  very  careful  reading  of  the 
original  manuscript  and  bis  helpful  comments. 
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